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We present the current status of our ongoing calculations of pseudoscalar meson decay constants for mesons that 
contain one light and one heavy quark (/s, /b 3 , }d, /d 3 )- We are currently generating new gauge configurations 
that include dynamical quarks and calculating the decay constants. In addition, we have several new results for 
the static approximation. Those results, as well as several refinements to the analysis, are new since Lattice '96. 
Our current (still preliminary) value for /s is 156 ± 11 ± 30 ± 14 MeV, where the first error is from statistical and 
fitting errors, the second error is an estimate of other systematic errors within the quenched approximation and 
the third error is an estimate of the quenching error. For the ratio /s s //s, we get 1.11 ± 0.02 ± 0.03 ± 0.07. 
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1. INTRODUCTION 

Heavy-light meson decay constants such as Jb 
and /b s are of great interest because they are 
necessary to interpret current and future mea- 
surements of B-B and B s -B s mixing. Only with 
knowledge of the decay constants and correspond- 
ing ^-parameters can we extract the Cabibbo- 
Kobayashi-Maskawa matrix elements and begin 
to decipher the fundamental information about 
the Standard Model (or physics beyond it) that 
is encoded in the CKM matrix. Recent reviews 
of this active field can be found in Ref . [Q . Here 
we will concentrate only on the recent work of the 
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MILC Collaboration. 

We have been involved in the calculation of 
pseudoscalar meson decay constants for some 
time now 0. Our initial goal was to calculate 
fs in the quenched approximation. This calcula- 
tion was done using gauge couplings between 5.7 
and 6.52. Our current goal is to estimate the er- 
ror that comes from the quenched approximation; 
however, we have also been refining our analysis. 

To estimate the quenching error, we have calcu- 
lated decay constants on lattices generated with 
two flavors of dynamical Kogut-Susskind quarks. 
The gauge coupling ranges from 5.445 to 5.7 and 
a variety of quark masses have been explored. 
(Table [l] summarizes the lattices that have been 
used in our calculations.) To go from an estimate 
of the quenching error to an actual calculation 
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Table 1 

Lattices used in our calculation. (Lattice Q is 
used only to study finite volume effects and is 
not included in the final results.) 



Quenched 



Name 


6/3 2 


size 


# 


At 


5.7 


8 3 x 48 


200 


B 


5.7 


16 3 x 48 


100 


Et 


5.85 


12 3 x 48 


100 


ct 


6.0 


16 3 x 48 


100 


Q 


6.0 


12 3 x 48 


235 


Dt 


6.3 


24 3 x 80 


100 


Ht 


6.52 


32 3 x 100 


60 


Dynamical, TV/ = 2 


Name 


6/g 2 , am 


size 


# 




5.7, 0.01 


16 3 x 32 


49 


G t* 


5.6, 0.01 


16 3 x 32 


200 


L 


5.445, 0.025 


16 3 x 48 


100 


N 


5.5, 0.1 


24 3 x 64 


100 





5.5, 0.05 


24 3 x 64 


100 


M 


5.5, 0.025 


20 3 x 64 


100 


P 


5.5, 0.0125 


20 3 x 64 


100 



t Approximately same physical volume 
** from Columbia group 
* from HEMCGC 



with dynamical quarks will require much work. It 
will be necessary to extrapolate in the dynamical 
quark mass to the u/d mass. It will also be neces- 
sary to extrapolate results to zero lattice spacing. 
The dynamical strange quark, which does not ap- 
pear in our calculations, will have to be included. 
Also, it would be nice to use the same lattice for- 
mulation of the dynamical and valence quarks. At 
the coupling 6/g 2 — 5.5, we have results for four 
quark masses. We plan to have the same number 
at 5.6, so we will be able to attempt extrapolation 
in the light quark mass for at least two couplings. 
We have no current plans to include the effects of 
a dynamical strange quark. 

There have been five significant changes be- 
tween the results we presented at Lattice 96 ||^] 
and those presented here. We now extrapolate 



to a physical n that gives the correct pion mass 
rather than to n c . We now evaluate the perturba- 
tive renormalization factors in a more consistent 
way. We have explored more fit ranges and this 
has helped us to refine our estimate of the statis- 
tical/systematic errors of the fits. We have new 
results for the decay constants in the static ap- 
proximation for several additional configuration 
sets. We are now able therefore to determine 
physical quantities at the 6-quark mass, on es- 
sentially all data sets, by interpolation between 
the static limit and propagating quarks that are 
lighter than b. (On two data sets, M and P, we 
still are forced to extrapolate at the present time; 
new static results on these sets will be completed 
shortly.) We have also refined our analysis of the 
7r mass and so that our chiral fits are better 
understood and controlled. Each of these issues 
is discussed in a subsequent subsection. We then 
present a summary of our results and, finally, our 
conclusions. 

2. CHANGES TO THE ANALYSIS 

2.1. Chiral extrapolations 

For the results presented at Lattice 96, all our 
chiral extrapolations were done to n c . We now 
determine the light quark mass from 

m n (Ki) = m^phys , 

where TO wp hy S is the physical pion mass, and / w 
is used to set the scale. This has resulted in 
a small reduction of fg of 2-6 MeV, for vari- 
ous cases. Physically, this change comes about 
because a/^ att is now bigger, which results in a 
larger value for a, and hence a smaller value for 
fg since afg changes less than af n . 

2.2. Perturbative corrections 

There are several improvements to the way 
we calculate the perturbative corrections to the 
renormalization constants of the axial-vector cur- 
rent. Tadpole improvement is applied in a con- 
sistent manner ||, and the remaining perturba- 
tive corrections are calculated using a boosted 
coupling constant g v (q*). The scale q* has been 
calculated for Z s ^ at by Hernandez and Hill to be 
2.18/a ||. For propagating quarks, the value 1/a 
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has been suggested by Lepage and Mackenzie [|| 
and this value is used in the analysis presented 
here. 

Table 2 

Differences (new — old) in values of Jb from cur- 
rent tadpole improvement 

Quenched 



6/g 2 difference (MeV) 



5.7 


+5-8 


5.85 


+3 


6.0 


+2 


6.3 


±1 


6.52 


-0.2 



Dynamical, Nf = 2 



6/g 2 , mass difference (MeV) 

5.445, 0.025 +3 

5.5, various 4-7 

5.6, 0.01 +2 

5.7, 0.01 -0.1 



Results presented at Lattice 96 followed an ap- 
proximate method of Bernard, Labrenz and Soni 
H) that was designed to reproduce the results of 
tadpole improvement at 6/g 2 — 6.3. However, 
that method differs from the more consistent cur- 
rent method, and the differences increase with the 
lattice spacing. The differences are shown in Ta- 
ble H We note however that a very recent lat- 
tice perturbation theory calculation by Bernard, 
Golterman and McNeile gives q* ~ 2.3/a [fjj for 
the Za for Wilson quarks. When this value is put 
into the analysis, the differences listed in Table |^ 
will decrease, making the analysis closer to that 
of Lattice 96. 

2.3. Alternative fit ranges 

To determine the pseudoscalar mass and decay 
constant, it is necessary to fit hadron propagators. 
Of course, it is not known a priori what range 
of distances to use in the fit. It is appropriate 
to consider the variation resulting from picking 
different ranges as part of the statistical error. At 
Lattice 96 we had a limited choice of alternative 



6/g 2 = 5.445, ani q =0.025 

0.6 — 




0.0 0.5 1.0 

1/M 



Figure 1. f M x M 1 ' 2 vs. 1/M for 6/g 2 = 5.445 
and am q — 0.025 (data set L) with no static re- 
sult. 



ranges used to estimate this error. We have had 
the opportunity to explore additional ranges, and 
this has changed some of our statistical errors. 
For example, for 6/g 2 = 6.0 we had thought that 
the variation from changing the range was ±15 
MeV. We currently estimate it to be ±8 MeV. 

2.4. New static fg calculations 

All the improvements described above can be 
implemented after the bulk of the computation 
is completed. However, the calculation of im- 
proved values for fs in the static approximation 
involved a reanalysis of many of our lattices. In 
our original calculations, the f^ atlc results were a 
byproduct of the propagating quark calculation. 
On large volumes, however, they were subject to 
contamination from states with non-zero momen- 
tum, and therefore were left out of the analysis. 
On small volumes the energy gap between the 
ground state and states with momentum 2n/L 
was sufficient to suppress the non-zero momen- 
tum states. 

Our recent calculations with dynamical quarks 
have been on larger volumes than the early ones, 
and we have done new calculations to determine 
the static decay constants. In some cases, we also 
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6/g =5.445, am q =0.025 




Figure 2. f M X M 1 / 2 ws . 1/M for 6/g 2 = 5.445 
and am q = 0.025 (data set L) including the static 
result. 

improved the accuracy of the quenched results. 
At Lattice 96, Craig McNeile [| described our 
method and some of the initial results. The key 
point is that the static quark is smeared relative 
to the propagating quark and the basis set in- 
cludes about 10 different smearing functions. The 
analysis of the static fs data is still being refined. 

With more accurate results in the static case, 
we are able to improve some of our fits for the 
decay constant vs. 1/M, where M is the heavy 
meson mass. Let's consider some examples. 

Data set L with 6/g 2 = 5.445 and am q = 0.025 
is a case where we had no previous static result. 
In Fig. we see fy/M plotted vs. 1/M. The 
curve is fit to the five decorated octagons. The 
burst is the extrapolated value of the decay con- 
stant and corresponds to 203(8) MeV. An alter- 
native fit to the points with 1/M > 0.5 gives 
200(7) MeV. In Fig. §, we see that including the 
new static value results in an interpolated value of 
210(7) for Jb- Using the points to the right, we 
find f B = 209(7) MeV. The results for the two 
choices of fitting range are more consistent and 
the error is slightly reduced. 

As another example, for case G with 6/g 2 = 5.6 




5 10 
distance from source 



and 



= 0.01 our old estimate of /J 



was 



265(12)(10), where the first error is statistical for 



Figure 3. Old static results on data set A for effec- 
tive masses and fitted masses, (a) smeared-local 
propagator, (b) smeared-smeared propagator. 



a fixed fit range and the second is the variation 
from the fit range. Our new result is 285(4) (1). 
The errors are considerably reduced and the re- 
sult is somewhat higher. 

Finally, we consider a quenched example, case 
A. In this particular case, although the volume is 
small enough to avoid contributions from higher- 
momentum states, we had poor plateaus in our 
effective mass plots. Figures || and [| show effec- 
tive mass plots for the old calculation and the new 
static calculations. In each case we show two ef- 
fective mass plots, one for the smeared source and 
local sink and one for smeared source and sink. 
The vertical bars near the bottom of the graph 
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Figure 5. Decay constant vs. the inverse meson 
mass on data set A. In the top graph (from data 
analysis presented at Lattice 96) the static result 
is ignored. In the middle graph the old static 
result is included in the fit. In the bottom graph 
our new static result is included. 



Figure 4. New static results on data set A for 
effective masses and fitted masses, (a) smeared- 
local propagator, (b) smeared-smeared propaga- 
tor. 



show the fitting range. The horizontal lines give 
the fitted mass. The confidence level of a fit in- 
cluding the old static result was just 0.09. With 
the new results it is 0.7. The new results ap- 
proach the asymptotic value more quickly and the 
errors in the effective masses are smaller. The old 
value for fg aUc was 216(27) MeV; the new value, 
283(18). Using the old static value to determine 
/b, we would have gotten 201(13) MeV. With- 
out the static value included in the fit, we ob- 
tained 227(18) MeV. (This was the fit chosen at 
the time of Lattice 96.) Including our new static 
value gives 229(12) MeV, in agreement with the 



Lattice 96 result, but with a smaller error. In 
Fig. ||, we see all three fits. The points determin- 
ing the curve are the decorated octagons. The 
bursts show both fs and fo- [Note: we quote 
statistical errors only for fs here. Choices of dif- 
ferent fitting ranges result in additional errors of 
order 15 MeV. However, the relations among the 
various treatments of the static point remain es- 
sentially unchanged.] 

2.5. Scale determination 

The final improvement we have made to our 
analysis involves the determination of the lattice 
spacing from /„. As mentioned at Lattice 96, we 
often had a difficult time getting good linear fits 
for m£ or vs. 1/k. We used quadratic fits to 
estimate the systematic error in the scale, result- 
ing in variations from 8 to 25%. However, our 
earlier fits for and f n were based upon a fixed 
fitting range for all three light k values. We have 
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Figure 6. Chiral extrapolation for and f n on 
data set D with a fixed fitting range 



relaxed this constraint and let the fitting range 
vary. This has certainly helped to improve the 
confidence level of the chiral extrapolation. As 
an example, for case D, we show both and 
/„■ in Figs. U and ^. The individual points in 
the plot are labeled with x 2 / (degrees of freedom) 
for the fits to the pion propagators as a function 
of distance from source. We see, particularly at 
the lightest k, that we are able to improve the 
fit by adjusting the range. The confidence level 
of the linear chiral fits is determined by the \ 2 
of the full covariance matrix. For the fixed fit 
ranges, we have confidence levels of 2 x 1 -7 and 
0.026 for linear fits to m 2 and f n , respectively. 
When we allow the fitting range to vary, the con- 
fidence levels improve to 0.20 and 0.56, for the 



corresponding cases. We have applied the same 
procedure to case O, and the linear fits in 1/k 
are improved, but not yet good. We have not yet 
implemented this for every case. 
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Figure 7. Chiral extrapolation for and / w on 
data set D with different fitting ranges 



We have both case A and B in order to study fi- 
nite volume effects on fs\ however, f n also varies. 
We have recently completed production running 
on a large scale finite volume study of the light 
quark spectrum at the same gauge coupling ||. 
We studied six quark masses on four lattice sizes 
from N s = 12 to 24. We should be able to study 
both the linearity of the chiral extrapolation and 
finite volume effects to higher accuracy. 
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Figure 8. fg vs. lattice spacing. 
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Figure 9. fB s /.fB vs - lattice spacing. 
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3. RESULTS 

Having carefully fit the hadron propagators and 
done the mass extrapolations, we are now ready 
to address the lattice spacing dependence of our 
results. We have results for fg (/d) and fg, 
Ud b ) and their ratios. Figure || shows our latest 
results for The diamonds are the quenched 
results and two fits are shown for their a depen- 
dence. The horizontal line is a constant fit to the 
three smallest lattice spacings. There is also a lin- 
ear fit to all of the diamonds. Both fits have quite 
reasonable confidence levels. The bursts plotted 
near a — are the extrapolated values. In Ta- 
ble ||, we list central values coming from the linear 
fit. The errors in the graphs are based upon sta- 
tistical and fitting errors, and are the first errors 
listed in the table. The second error is our esti- 
mate of the systematic errors within the quenched 
approximation. The chief error here is that from 
the ambiguity in how to extrapolate in a. Other 
sources of error include the chiral extrapolation, 
finite volume effects, weak coupling perturbation 
theory, the extrapolation in \/M and large ma ef- 
fects. How we deal with each of these effects has 
been discussed in Ref. (2). We must admit that 
we are wavering between taking our central value 
from the constant and the linear fit to the lattice 
spacing dependence. Since we are currently us- 



ing the value from the linear fit, the second error 
is more likely to be positive than negative. The 
error from using the quenched approximation we 
consider to be rough at present. We now have 
results for seven sets of dynamical quark param- 
eters. They are denoted by crosses in the figures. 
When we have finished our calculations with cou- 
pling 5.6, we will attempt an extrapolation in the 
dynamical quark mass. 

In Fig. |[ we show the ratio JbJ Jb as a func- 
tion of lattice spacing. One expects that in a ratio 
many of the systematic errors will be common to 
the two quantities and that the ratio can be calcu- 
lated more accurately than either quantity. Here 
there is better agreement between the constant fit 
to the three weakest couplings and the linear fit 
to all the quenched couplings. There is one point 
that seems to stand out from the others. This is 
case G, an analysis of I6 3 x 32, 6/g 2 = 5.7 dy- 
namical quark lattices that were provided to us by 
the Columbia group . These lattices are phys- 
ically the smallest lattices we have studied. We 
tried to test whether we might be seeing a finite 
size effect by doing a new quenched calculation on 
a small volume. Our most comparable quenched 
coupling is 6.0 where we studied 16 3 x 48 lat- 
tices, so we reduced the spatial size to 12 for this 
calculation (lattice Q). This smaller volume cal- 
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Table 3 

Preliminary results for meson decay constants 



fB 


= 156 ±11 ±30 ±14 MeV 


fB. 


= 173 ±9 ±40 ±18 MeV 


fsjfs 


= 1.11 ±0.02 ±0.03 ±0.07 


fD 


= 192 ±9 ±18 ±10 MeV 




= 209 ± 7 ±27 ±12 MeV 


Jd 3 1 Id 


= 1.09 ±0.02 ±0.05 ±0.05 



culation is plotted in Figs. |8j and |9j with the fancy 
diamond. Although it appears a little lower than 
the larger quenched volume, the effect is not as 
pronounced as for the Columbia lattices. Thus, 
we do not have a credible explanation of this as 
a systematic effect. 

Of particular note among our results is the 
value of fu s . The world average experimental 
result has recently been reported in a review talk 
by J. Richman 0. The result, 241 ± 21 ± 30 
MeV compares quite reasonably with ours. 

4. CONCLUSIONS 

We feel that our calculations of the meson de- 
cay constants have progressed to the stage that we 
can control the sources of systematic error within 
the quenched approximation. At this point, the 
major source of uncertainty within that approx- 
imation is whether to use a linear extrapolation 
for all of our couplings or a constant fit to the 
weakest ones, to extract the continuum results. 

Our calculations with dynamical quarks are 
now giving us a rough bound on the error from 
quenching, but we are not yet in a position to 
extrapolate in dynamical quark mass or lattice 
spacing. We are still running on lattices with 
6/g 2 — 5.6. Beyond these calculations, one can 
imagine putting in a dynamical strange quark, 
or repeating this calculation with dynamical Wil- 
son quarks. Such ambitious plans are not on our 
agenda at the present time. 

The apparent agreement between lattice calcu- 
lations of fu a and the recent experimental results 



is encouraging, but the errors on both arc still 
rather large. With new experimental facilities in 
Japan and the United States coming on-line in the 
near future, and with continuing improvement in 
the lattice results, we expect increasing interplay 
between experimentalists and the lattice commu- 
nity. 
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